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Twin beam fluctuations are analyzed for detuned and mismatched OPO configurations. Reso- 
nances and frequency responses to the quantum noise sources (quantum and pump amplitude/phase 
fluctuations) are examined as functions of cavity decay rates, excitation parameter and detuning. 
The dependence of self- and mutual correlations of beam amplitudes and phases on detuning, mis- 
- match and damping parameters is discussed. 
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1—5 . I. INTRODUCTION 

in : 

Continuous wave twin beams owe their popularity to the mutual quantum correlations. It has been shown theoret- 
ically (see f.i. [0-|8)) and proved experimentally (see f.i. |^|,|o|) that this feature can be exploited for replacing classical 
fields in shot noise limited optical measurements. In particular, the spectrum diffcrcncc of the two beams is shaped 
like a Lorentzian with a rather good noise suppression at zero frequency. On the othcr hand, Fabre et al. ]ïï| ] havc 
called the attention on the squeezing exhibited by each beam for pump intensity larger than four times the threshold 
■ value. They have also shown that additional squeezing can be imposed on each beam by using the photocurrent 
ÇSJ \ detected with the other one. 

These properties are contrasted by unsuited damping coefficients of the OPO cavity, imperfect tuning of the three 
mode resonances, limited intensity and amplitude/phase fluctuations of the làser pump. First Lane et al. p2[ drew 
the attention on the critical dependence of the noise suppression in the difference spectrum on the cavity mismatch, 
as experimentally confirmed by several authors (see f.i. fL^JÏ4|| ). Stable and narrow line làsers emit in general low 
Oh, intensity beams. Then, for observing single beam squeezing it is essential to lower the threshold by reducing the pump 
+1» ■ mode decay rate. This solution makes the OPO a triply resonant device featuring relaxation oscillations Jl5|— p~T|] . Due 
to the unavoidable OPO crystal losses the amplitude/phase fluctuations of the pump contribute to the difference 

d ' 
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spectral noise at low frequency, where quantum noise suppression is expected. The situation becomes worse and 
worse as the cavity is detuned from the resonant configuration. 

These few remarks point to the several parameters which the features of an OPO depend on, namely, cavity damping 
coefficients, mismatch of these parameters, degree of excitation above threshold, crystal absorption, deviation from 
k> , resonance condition (detuning), and pump amplitudc/phasc fluctuations. The last ones play a prominent role in 
5_j ■ the lower part of the frequency spectra of all relevant quantities. Many papers have been published discussing the 
| role of some of these parameters. Here we aim to discuss analytically the twin beam fluctuations by encompassing 
systematically the many parameters mentioned above. 

The representation of the three modes (signal/idler k = 1, 2, pump k = 0) in the form au — r^e~ % ^ k (I + /ifc) 
lends itself to separate the steady state amplitude Tk from the fluctuating relative amplitude /i^ and phase <pk- In 
regard to the fluctuations, the OPO behaves as a forced linear system. Then, the Fourier transforms àk (w) of the 
combination 5a,k = [ik — i(f>k and the adjoint à\ (u>) = à\ (— u>*) form an algebraic system. The dependence of the 

àk (oj) on the noise sources is embodied in the frequency responses KJ, 2 (u>) of the k-th mode to the Z-th noise source. 
These responses are analytic functions of u) with, all together, six complex poles (system resonances) in the half-plane 
Im (u>) > 0. Their location (gcncrally complex or purely imaginary) influences the field correlations in the time domain 
and, hence, the relative frequency spectra. For example, when the poles have large real parts the system exhibits 
relaxation oscillations, responsible of typical peaks in photocurrent spectra. In short the K k l 2 (lo) and their poles 
provide a link among the OPO, pump làser parameters and the statistical features of the twin beams. 
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While several authors have shown that the unavoidable cavity mismatch reduces the mode correlation, less attention 

has been paid to the effects of detuning, apart from the analysis of Ref. Jïï| dedicated to OPO with pump bandwidths 

much larger than the signal/ idler ones, i.e. in the adiabatic limit. For large detunings bistability, self-pulsing 

and chaos may appear |l§|Jl9|] , while for moderate deviations from cavity resonance, modifications of the spectra of 

the single beams and their differences are expected. In presence of detuning the decay constants become complex, 

K-k = Jk + i<Pki and the equations of motion of /ik do not separate from those of 4>k- In particular, the pump phase 

diffusion leaks into the single mode amplitude spectrum by enhancing the disturbing effects of the classical pump 

fluctuations. From the analytic point of view, the detuning has the effect of breaking the symmetry of the frequeney 
o — o — * 

responses, K,) 2 (ui) ^ KJ, 2 (— ui), so that /ifc (t) and 4>k (í) do not commute with the respective quantities at different 
times. 

The model, herein discussed, is an extension of the quantum analysis of Graham and Haken [^,^ÏJ and Fabre et 
al. pH| . In previous works the variance and covariance of the OPO beams have been analyzed by representing the 
fields either as operators jÏ5| or classical variables. While Reid and Drummond |2^] have taken in full account the 
fluctuating phase factors e~ l ^ k by using the stochastic equation method, the operator representation has been limited 
to models neglecting the phase diffusion processes. To fill the gap between these two approaches, quantum Langevin 
equations will be used for calculating second and fourth order correlations. The operator representation offers the 
advantage of introducing the commutators of the field variables at different times. Some of these commutators vanish 
at resonance, so that their Fourier transforms provide an additional signature of the detuning. Essential to these 
calculations will be the assumption of Gaussian statistics of the mode operators. Most of the numcrical examples 
discussed below refer to OPO exhibiting relaxation oscillations, that is having pump decay rates comparable with the 
signal/idler ones. 

The paper is organized as follows. In Sec. 2 the bàsic equations for an OPO are reviewed. The steady state 
solutions are recalled in Sec. 3, while in Sec. 4 it is introduced the linearized set of the equations of motion for the 
fluctuating parts of the infernal modes. Then, switching in Sec. 5 to the frequeney domain, the Fourier transforms 
of the fluctuations are expressed as functions of the noise inputs. Moreover, the quadratures of the single beams are 
represented by combinations of the quantum and classical (pump amplitude and phase) noise sources times suitable 
transfer functions. The system resonances are discussed in Sec. 6 by dwelling on the dependence of the pole location on 
the cavity decay constants and excitation parameter. Eventually, the photocurrent fluctuation spectra are examined 
in Sec. 8 with the help of the transfer functions K,l 2 discussed in Sec. 7. 

Details on the solution of the algebraic system describing the fluctuations of the OPO modes in the Fourier domain 
are given in Appendix A and B. The damping coeficients of the OPO cavity and the losses of the nonlinear crystal 
have been indicated respectively by jk and k^, by adding an apex for the total losses 7Í = 7fe + >ffe- Quantities 
normalized with respect to the average damping 7 = (7J + 72) /2 have been represented by superposing a tilde. 
<5i = (7J — 7 2 ) /2 = —82 measures the cavity mismatch. 



II. EQUATIONS OF MOTION 



The evolution of the mode amplitudes üq (pump), a± and a 2 (signal/idler) of an OPO is described by a system of 
Langevin equations 

q í ai \ / -K[ai + 2xa 4 + R'i \ 

-Q t £12 = - K ' 2 a 2 + 2 X a Q a\ + R' 2 (1) 
\ a o / \ -K , a -2x*a 1 a 2 + e + R' ) 

with R' k (t) independent fluctuating delta-correlated Langevin forces and e = |e| e _i ^ p with the pump phase <j) p (t) 
undergoing a slow diffusion, ^(</> p (t) — <p p (O)) 2 ^) = 2Avlt. The coeficients n' k (fc = 0, 1, 2) stand for 

K k = lk + Xk+ Wk = 7fc + Wk = Kk + Xk (2) 

with 7^ and >Ck the cavity mode damping rates, associated respectively with the output mirror and the other loss 
mechanisms (mainly duc to the non-lincar crystal), ip^ the detunings and 2x the strength of the parametric interaction. 

The operators R' k are proportional to the vacuum fluctuation òm fc ,i entering the system through the output mirror 
and to the noise &m fc ,2 associated with the crystal losses, 

R k = \/ r 2^fkbin kl l + V2xfc6m fc ,2 = Rk + V2Xkkn k ,2 (3) 
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Applying now the time reversal to the above Langevin equations ||,^2|] it can be easily shown that output, internal 
and input fields (b ou t k .i, a-k, &m fe .i) & r e connected in presence of detuning by the Standard boundary conditions fï~ïf| 

bout k ,i = \Z^ïka k - b inkA (4) 
Finally, the amplitude-phase representation 

a k = r k e-^ (1 + fik) (5) 

will be adopted in the following, with r k the steady state amplitude, <p k the phase, and [i k the relative amplitudc 
fluctuation. 



III. STEADY-STATE SOLUTION 



Dropping the Langevin forces and the fluctuating amplitudes, the system (Q) admits the steady state solution (see 
Eq. (g) 



a k = (a k ) = r k e l<tn > 



(a* k = Although (p k fluctuates it has been included in order to get average vàlues of some phase combinations. 

The signal/idler amplitudes aj (j — 1,2) are different from zero for |e| above a thrcshold 



6 ~ 2 X ~ 2 X cos^ - |Ko|r ° (6) 

only if the relative detunings tpj (see Eq. (^)) are proportional to the damping factors 7^ = 7^+ >cj of the respective 
modes, i.e. 

fi <P2 . , <PO.fi , s 
— - = — - = tan tp , — = tan^o - V) 

Ti 72 7o 

Most OPO cavities are stabilized on the pump mode, so that in the following we will put ïp = ipQ. 
Amplitudes r k and phases (p k of the three modes result interconnected 

Wx\tI = W 2 \tI = {£-1)\k' q \tI 

(01 +02-00) = ïp (7) 

with £ related to the amplitude parameter E — \e/e th \ by 



£ = \J E 2 - sin 2 ïp + 1 - cos^ (8) 

In absence of detuning £ reduces to E. Notice that the relation between the frozen phases of the three beams (see 
f.i. Eq. (2.9) of Ref. p^fl ) depends on the detuning through the phase ïp. 
In turn tq and (po are respectively given by 



r =~ 



2X 

60) = j> ~ Í> v (9) 



with <p p the phase of the external pump and tp p = arcsin (smip/E). 
In the following the parameter 

C 2 = (£-l)W \rl 

will be uscd. 
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IV. FLUCTUATIONS 



The phase relations obtained for thc steady state solutions are satisfied only approximately. In particular, the phases 
!>i + 4>2 — 4>o and 4>o — «fip (Eqs. (0-b) and (§)-b)) undergo stable fluctuations, while <fii — </>2 follows an undamped 



diffusion process. On the othcr hand, 



a k = r k e 



°Sà k 



where Sa k = pik — i<fik (Eq. 

By exploiting the smallness of the quantities <pi + <p2 — 4"3 ~ V' and 
linearized with respect to the fluctuating phases 4> k and amplitudes \x k 



(10) 

- (f> p + ip — ip p the system (|j) can be 
Next, changing the phases (j>k and 4> p 



mto 



^i/) p , (f>' = (f>o + ip — ip p , and 4>' p = (f>p + tpp does not change the statistical properties of the system 



solutions. Thus, keeping the old symbols for the new quantities, .the linearized system reads: 

Sài 



- — h 5a\ — Sü2 — Sao = Z[ 
n i 

—y- + (5ct2 — Sao — Sai = Z' 2 



Sào 



- — - — - + 5a + Sai + Sa 2 = Z' 
£ - 1 V «o / 



with 



Z'o 



Ze 
Za, 



(i = 1,2) 



1 



£ - 1 
E 



Koro 



£ - 1 
E 



£ - 1 



having denoted by /i e = Se/e and p the pump relative excess noise and phase fluctuation respectively. 
The output field, having average intensity Ij = 27-,-r 2 = (2jj/j'^ cosV> C 2 , is given by (see Eq.([|)) 



(11) 



(12) 



(13) 



where Zj is given by (|Ï2|-a) with R'- replaced by Rj . 

Before concluding this section, it is worth looking at the approximations underlying the linear syste m (0 ). Products 
of phase-factors e l ^e" lx can be rcprcsentcd by thc gcncralized Campbell-Hausdorff formula (see f.i. p4| p. 35) 



e ** e -*x = exp Utf, - i X + - [4>, x] + j2 [[<!>> x] » 4> + x] ■ ■ ■ 

so that choosing = 4>(t), x — 4> (t + dt) we have for the derivative of a k (t) 



à k (í) = r k e- 1 ^ ( Sà k (í) - - \4> k (t) , fc (í) 



Vk,<Pk 



which reduces to Eq. (HO) at the first order in the fluctuating quantities. 



V. FREQUENCY ANALYSIS 

Now it is worth introducing the Fourier transforms 
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sa k m i rwM^)^ (14) 



5a[ (t)J V2tt J-oo V«fc M 

with à\ (ui) = à k (— u>*) the adjoint of the Fourier component relative to the frequency — ui. In the following we will 
consider generally complex freqüències u> and the symbol |, operating on the product of a c-function times an operator, 
will indicate the transformation £> (ui) à (ui)] = /* (—ui*) (—w*). 
In the frequency domain Eq. (1ÏJ) reduces to the algebraic system: 

Aiòi — Úq — à 2 = Z[ 

A 2 à 2 - ào - à\ = Z' 2 (15) 

A à + àx+à 2 — Z' + Z f + Z^ 

with 

lj f ui\ 
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According to Eq. (12-c,d) Z í — -^-^e l ^ v \x t and Z$ = — i-^\ eV ^ p 4>p depend respectively on the pump amplitude excess 

noise pL t and on the phase cf> p . 

Next, completing the system ( |Ï5|) with the respective |-transformed equations we obtain for the Fourier transforms 
àh (ui), and à\ (u>) a sixth-order algebraic system depending on the forcing terms Z' k , Z e , Z^, Z k , Z\, z|. Being 

an OPO a stable system (see f.i. Ref. p^|), à k (ui), à\ (ui) may become singular only in the upper imaginary plane 
(lm(w)>0). 



A ^ = 1+ ^< Ao = 1+ï z —i (16) 



Now, solving the system (|15| ) completed by the J-counterpart yields 

\\C k - B k C\ 



«fe 



A k A{ - B k B k 



with the functions A k (ui) , B k (ui) , C k (ui) defincd in Appcndix A. 

Next, introducing the Fourier transformed quadratures of the output modes p,p h ~ | (jí k + pj^j , and 

§ (jík — (ï]^j we have (see Appendix A) 

where A ® A" is a shorthand for | ( A (w) 6 (w) + A* (— w*) 6* (w)) . 

The functions K k ] 2 represent the frequency responses to the quantum noise sources and the classical phase and 
amplitude fluctuations of the pump lascr. Thcy arc analytic in ui with poles in the upper half plane Im (ui) > 
coincident with some zeros of AfcAji — B k , as shown in the following sections. 

According to flÏ7j ) the amplitudes and phases depend on the quadratures X® — | [e %e Z k + e~ lB Z^j of the fc-th 

vacuum noise entering the cavity through the coupling mirror and that relative to the crystal loss x k , X x . The 
fluctuating amplitude /Lt e and phase <fi p of the external pump are represented respectively by Xf and A^ . 



(17) 



VI. RESONANCES 



Important features of the OPO fluctuations can be investigated by extcnding the system (15) to complex vàlues of 
ui. This will be done by retaining the above definition of the operator J transforming an analytic function of u>, with 
poles Q r , in a new one with poles — Í2* and complex conjugate residues. The system of six unknown functions à k (ui), 



Sí (u>) is characterized by the sixth order characteristic polynomial D (ü>) = D' (oj) ü with a zero at the origin, 
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D' («) = Q 5 + D^Cj 4 + D^Cu 3 + D^ú 2 + D^ú + D« 



(18) 



with ió the frequency normalized to 7. The coeficients are reported in Appendix B (Eq.(|37|)). 

For a detuned balanced cavity (5 = 1 5x,2 1 = 0) Z)W are independent of tan-0 so that D' coincides with that of a 
resonant OPO having an effective excitation parameter E e f f — 1 + (£ — 1) / cosip. In particular, D' factorizes into the 
product D + D_ of polynomials discussed in the following for the resonant case. This means that a detuned balanced 
OPO cannot be distinguished from a perfectly resonant one with an appropriate pump strength. 

The six zeros fl r (r = 0, . . . , 5) of D (w) represent the system resonances. In particular, the real parts of their vàlues 
indicate the oscillation freqüències, whereas the imaginary ones stand for the corresponding damping parameters, i.e. 
the resonance widths. Since D (üi) = D* (cD), each resonance fl r having a finite real part is mirrored by — fi*. Apart 

from the one at the origin (fl = 0) , the resonance ííi = iwi is located on the positive imaginary axis, while the 
remaining four ones may be either purely imaginary or may have finite real parts. In fact, for excitation parameter 
varying in an interval depending on the cavity damping factors, one or both couples fl r , —fi* are replaced by two 
purely imaginary roots. In the case of non-vanishing real part fl2 = — fl% = ^2 + ÍW2 and fl^ = — fit, = £4 + ÍW4; the 
system is characterized by two dampings u>2, 104, and freqüències Q2, In particular, w\ + 2w2 + 2u>4 is independent 
of the detuning, whereas for a balanced OPO w\ — 2. 

The expressions of fl r are rathcr lcngthy and complicated, thus, for the sake of simplicity, we briefly discuss their 
major features, and dwell on two noticeable cases. Fig. 1 displays the dependence of the real and imaginary parts of 
fl r on the pump parameter E for ip — 0.3, 7Ó = 2 and a 20% mismatch (5 = 0.1). The resonance freqüències ü>2 and 
ÒJ4 start from the origin and grow approximately like \fE. The damping parameter w\ starts from zero and rapidly 
reaches the constant value of 2, whereas W2 = 2 is independent of E 7 and u>4 moves quickly away from 2 toward 1. 
This behavior is not very sensitive to the detuning angle. 

At resonance D' factorizes into D+D-, being D_, and D + two polynomials of second and third order respec- 
tively, 

L>_ (w) = Cò 2 - i 7 (1) w - 2 7 (2) , D + (u) = w 3 - i 7 (1) w 2 - 7 (2) w + i 7 (3) 

where 7 « = 2 + %, 7^ = 2E%, and 7^ = 4 (E - 1) % (l - 6 2 ^ . 

D- has two roots 04,^5 which are purely imaginary if 7W 2 > 87^, otherwise O4 = — fit,- In turn D + has three 
purely imaginary roots (Oi = iw\, fÍ2 = ÍW2, and ÍÍ3 = ÍW3) if 

7 (1)2 > 3 7 (2) j 7 (1) 7 (2) > g 7 (3) ) ^(1)2 47 (2)^j 7 (2)2 > 27^)2 + 3 Aj^ 1 ) 3 - CvyM^ 2 )) 7(8) 

According to these inequalities (derived by means of the Sturm sequence) the roots are all purely imaginary for 
1 < E < E max , with -E max depending on % and 5. In the example plotted in Fig. 2 for 5 — 0.05, E max grows with 
7q except for the presence of a peninsula in the interval % ~ 1.5 -r- 2.75. Scaling the threshold with % (Eq. (|h), 
for 7g ^> 1 the roots are purely imaginary for all realistic vàlues of the excitation parameters. This explains why the 
relaxation oscillations (resonances with non-vanishing real parts) can be observed only for 7 not much greater than 

7- 

While the poles of K? k control the transient behavior of i^p' 2 (r, 0)), the relaxation oscillations correspond to those 
of K® k . In particular, for E rather large 

Òi ~ i (l - 5 2 ) fl 2 ~^(l + ó 2 +%)+^2E% (19) 

Similar expressions were obtained in Ref p6| on the basis of a rate equation model for the relaxation oscillations. 

In Fig. 3 we have plotted the real and imaginary parts of ^2,3 versus 1 < E < 10 for different decay rates 
% = 0.5, 2.5, 4.5, and 6.5 and S = 0. These plots also display the dependence of the roots on the detuning angle ip of 
a balanced system having an effective excitation parameter E e ff equal to E. 

In the adiabatic limit (7Ó 3> 7 , see f.i. Ref. |ÏI]]) D' (ó>) reduces to a third order polynomial, 

D' (w) = ü 3 + D^Co 2 + D^Cj + D™ 

that is the OPO is characterized only by four poles, fio = 0, fl\ = iw\, fÍ2 = —fit, = fl2 + ÍW2 and the relaxation 
oscillations disappear. 
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VII. TRANSFER FUNCTIONS 

The K° k \ 2 are singular in the half-planc Im (çij > 0, their polcs coinciding with some rcsonances Ú r of the whole 

system. Their expressions for the general case, together with those for the adiabatic case and for a resonant device 
are reported in Appcndix A. 

For M ->«>,*£*-> - K 7f - X^T and < f Ak- 

3 3 3 ' 

K® k is not singular at the origin contrariwise to K?- and K?j, . In particular in proximity of the origin 



1 / tan^ 
1 i- 



KJ. ( W ) « ( W ) 



. 1 - ó 2 1 



4 cos 2 ?/> üj 



f (£ e// - 1) " (tan 2 t/» - (ggAf - 1) + ige/j tan V>) Sj 
K j0 M — (20) 

where (E e ff — 1) = {£ — 1) / cosí/; represents the effective distance from the threshold of a detuned device. 
Analogously the phase reduces for wsíOto 



1-à 2 1 /^, ~, \ (ge// - 1) - (tan 2 ?/> - (£ e// - 1) + iE eff tan V>) ój ~ f 



(21) 



namely the detuning enhances, through the factors cos -0, the phase fluctuations associated to the quantum noise 
sources. 

The mathematical difficulties associated with the singularity of K?j , K?j, at the origin can be avoided by displacing 
slightly the zero íl = t-^> íl = ís (e > 0), and letting e — > once performed the integrations. 

The presence of poles with real parts larger than the imaginary ones is evidenced in the plots of I I , l-K'?^! and 
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\KjQ\ reported in Fig. 4 for cavity parameters 7q = 1, S = 0.05, and E = 2, 3, 4. They exhibit a behavior typical of 
triply resonant devices, namely peaked at freqüències increasing with the pump rate. While the peak height of the 
self-coupled response ÜT?- (Figs 4a) slightly decreases with E, those of the cross-coupling K^-, (Figs 4b) and of the 
response K® (Figs 4c) to pump fluctuations keep, respectively, nearly constant and growing in height with E. All 
thcsc frcqucncy rcsponses are not very sensitive to the detuning. 



VIII. SPECTRAL DENSITIES 

For vacuum input fields, only antinormally ordered terms contribute to the variances, i.e. 

Conscqucntly, the noise corrclations have in the frcqucncy domain the following expressions: 

1 \ K j\ 



z e h z\ = (z* M z t m)' = s £ (u) 

( W ) Z\ (u/))' = (Z\ (w) Z M)' = (^^j S (w) (22) 



all the other terms vanishing identically. S t (o>) = (J£i € (u>)\ j and (uj) = Av^/u) stand for the spectral densities 

of the pump amplitude and phase respectively. The apex (()') indicates the omission on the right-side of the factor 
S(lü* +lü')/C 2 . 
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A. Gaussian statistics 



As a result of the superposition of uncorrelated spectral components, the single beam fluctuation amplitude Sa (t) = 
ka is proportional through a coefficient k to the annihilation operator à of an oscillator associated to a Gaussian density 
operator (see p5| , p6[ ) 

p = J^ e ~(n fí +n 4 ,)K + ""l n · t ' (K++K-) _ J^f^uKo £ vK + e ~wK- ^23) 

with 2A + — a* 2 , 2A_ = ò 2 , 4A'o = à'a + àà' forming a realization of the SU(1,1) group, with the relevant 
algebraic structure displayed by the commutators [A_, A + ] = 2Kq , [Kq,K±] = ±K±. In the above equation p has 
been expressed either as a single exponential depending on the coemcients n M and or as a disentangled product of 
exponentials depending on u, v and w. The connection between these two representations is displayed by the relations 

sinh .JUané sinh (■Jn ih n6 + 6) sinh Jn^n^, 

2u = In (—w/v) , v = q , w - 



sinh 2 6 sinh (,yn M n0 + 9) 

with coshO = (rifi + n^) I (n^ — n$). Here we limit ourselves to report the results of calculations based on the 
Wei-Norman technique (27| (see also |2^] ) which will be discussed elsewhere together with the calculation of the trace 

Tr (e 2uKo e vK+ e- wK ~) 



y/(vw + 1) e u - 2 + e~ u 

The moments of Kg and K± can be easily obtained by diffcrentiating the operator e 2uK o e vK +e -wK- w jth respect 
to u, v and w. In particular, 

|[«= i fM4(ín)·w4(í + ï) 

It is straightforward to extend the above representation to the twin pair 8a\ (í) = kà, 5ci2 (í) = kb for a balanced 
OPO, 

p = AAexp (C +A+ + C -iC - + C+ (K+ + K+) + C~ (K+ + Al)) 

with 4if^ = {àVà} 2 + i 2ií í = A 'oa + K ob , 2K± = {a&/a+&t j an d 4A~ = à f ò + àò f forming an enveloping 

àlgebra of SU(1,1) characterized by the algebraic structure 2 A,^, = ±K± , 2 [A:* 1 , A^:] = KZ all the other 

commutators vanishing identically. The disentangling of this operator will be discussed elsewhere. 

Using the squeezing operator S (6) — exp [| (à 2 — <^ 2 )] with tanh (26) = 1/coshO, introduced by Gardincr and 
Zoller p5| for transforming a gaussian density matrix into a thermal oscillator, we obtain in the present case p — 
2 sinh v "m"-» ex p (_2^/rï M rí0Ao) . Accordingly the twin beam fluctuations can be represented as a pair of squeezed 
and coupled thermal oscillators. 

Finally, from the Gaussian statistics it descends for the phase factors e 1 ^ 1 ") correlation 

i<t>j(-r) p -i^(0)\ _ „<(0 J (r)-0 J (O))0 3 (O)) 



B. Correlations and frequency spectra 

The various correlations can be represented as combinations of the functions 



K^4n (r) = E & E ^^M* (K) 

k r=l 
5 

M = E & E e in ^A-° fe , r < (0?) 
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w = E ^ E énrTK %,rKf k * ík) 

k r=l 
5 

(r) = e E e<nrr ^l,r^ 



(25) 



with Z = 0, 1, 2 , the sums being extended to the 5 roots Í2 r of the characteristic polynomial D' (i.e. except for ÍIq)> 



and KA" r standing for the residues of Kj' k 2 (tu) at Q r . 
We have for the correlations 

((&(r)-&(0))&(0)) 
( Mj (r) W (0)) 
( W (r) 0, (0))-<^(r) W (0)) 



4C 2 

4C 2 
1 

4C 2 



(Km, (t) - (0)) - Ai/j 

K NN (t) 

(r)+K^ m (r)) 



(26) 



and the commutators 



([0i(O),^(r)]) =^Re{^,(r)} 



([ W (0) lW (r)]) 
[0í(O),Mí(t)] 

([<Mo),Mt)]) 



2C 2 

-^ Re {^ 3 (r)} 
— Im{^. M . (t)} 



^ 

2C7 2 



(27) 



with t > 0. While (|27j-a,b,c) vanish at resonance, the last expression does not. For a detuned OPO fij (t) and <f>j (t) 
do not commutc with the respective quantities at different times while fij (i) does not anticommute with <j>j (t'). By 

similar reasoning it can be also shown (^(f>j (0) <f>j (0)^ ^ 0, that is the instantancous frequency is correlated with the 

phase. In Fig. 5 the phase commutator has been plotted versus r for an OPO with 7q = 4, 6 = 0.2 and ip — 0.1 and 
£ = 1.5, 4 (Eq.(||)). These plots exhibit respectively oscillatory and damped behaviors since for £ = 4 the roots have 
non vanishing real parts, while in the other case they are almost purely imaginary. 
The variance of the phase delay is given by 



with 



(r) - 0, (0)) 2 ) = ^Im {Km (r) - (0)} + 2Au jT 



Afj = Fa,Av l 



1-5' 



8C 2 cos 4 ip 

standing for the broadening factor of the pump linewidth Az/£, (see below Eq. (|ï|)), 

(E e ff - 1) - (tan 2 ?/; - (E e ff - 1) +iE e fftaaip) 5j Ecosip 



F ó = 



E, 



eff 



ZEeff 

In a similar fashion we have for the phase sum variance 

((& (r) + cf> ó (0)) 2 ) = ^ Im W " K%.,^ (r) + (0) - K^., (0)} 

so that (see [|§) in view of Eq. (|3) 



2Ai/jT 



,í^'(t-) p -%(0)^ 



- Al/, r 



+ 4^2 (r)-A'o,^(0))) 



e-^^e-^'W) = e " A ^ 1 + ^ (i^., (r) - tf^, (0) 



9 



while the second-order correlations for the mode amplitudes a,j , üj read 



( O t(r) Oi (0)) =r)e~^ 



1 + ( K HH W + W - ^ (0) - K Mj<l·l (r) - (r)) 



+ ^2 ( T ) + ( T ) " K Wi' (°) - ( T ) - W) 



Further, the intensity variances of the single beams are represented by 



a) (r) aj (r) , a) (0) (0) :) = r| (2 ( Mj (r) ^ (0) + Mj (0) Mj (r)> + i ([^ (0) , ^ (r)]>) 



,t 



(28) 



standing for time-normal ordering. Analogously, for the respective spectra, normalized with respect to r^/2C 2 , 



Sj + S 



with Sj separating in three terms 



E 



£-1 



E 



£ - 1 



Sé Sé 



(29) 



representing in the order the quantum and the pump amplitudc and phase contributions 

S N (w) =a% M + o^t-uO 

Sfa (w) = a\ } (ui) + al. (-lü) 



(30) 



S e and ^ represent the pump amplitude and phase spectral densities (see Eqs.(p2|-c, d)) while 0^ ,0^. and cr? stand 



for 



^-2 \ " I j^O I 2 J 



pitpp jy-0 _ p -v4> p TsOt 



In turn, the normalized spectrum of the commutator <pj] splits in the components 

S M]j (W) = ^,0] . M + <7p Mi0] . (-W) 

with 



(31) 



(32) 



Similar expressions hold for the output beams with Ij = 27 J r 2 in place of r| and cr^. (Eq. (|3T|— a) ) replaced by 



^2 \ " I 7^0 I 2 I 



K 



K 3 

jj 



27, 



I mI 7 / IS 



(33) 



In Fig. 6, we have plotted «Sj (Eq. (p9|)), the spectrum measured by analyzing the photocurrent with a spectrum 
analyzer, for different detunings (ip — 0,0.25,0.5) and pump parameter E = 3 by omitting the làser excess noise and 
phase fluctuations. For both 7Ó = 4 (top) and 0.5 (bottom), in the case of zero and small detuning, Sj starts from 
below the SNL, while, for tp = 0.5, it starts just above the SNL. In all cases Sj rapidly increases by reaching a peak in 
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correspondence of the relaxation oscillations frequency, and then decays monotonically toward the SNL. The height 
of the peak increases notably either with the detuning or by reducing the pump mode bandwidth. 

Figure 7 describes the effects on the single beam spectrum of the pump excess noise at different detuning angles. 

Íe~~ï) Sej ( w ) nas been plotted for different detunings (-0 = 0,0.25,0.5), 5 — 0.05, pump parameter E = 3, and 
crystal losses different for the signal/idler and the pump 2*1,2/7 = 0.3, xq = xi/3. For both cavity losses % = 4 
(top) and 0.5 (bottom ), (^^zr^j S ej increases with ip and u> by reaching a peak in correspondence of the relaxation 

oscillation. The behavior observed for 0.5 < % < 4 leads us to conclude that for mode dampings of the same order, 
the detuning confines the observable squeezing to a small fraction of the twin-beam bandwidth. Moreover, it increases 
the amplitude of the excess làser noise transferred to the single beam spectrum, thus preventing the reduction below 
the shot noise level. 

The difference spectrum Sd is obtained from Sj by replacing in (|3Ï|), ( p2| ) and ( |33| ) K- k 2 with AK k ' 2 = 
ïiK® k 2 — L·K^ 2 J I {Jx + Ï2) and analogously for the residues, ^f- and 2f. In particular, the contributions of 



pump amplitude and phase fluctuations are represented by Eqs.(|30|-b,c) respectively, with a 2 . and cr^. replaced 
by a 2 d and cr^ d , which are in turn defined by equations similar to (|3Ï|~b, c) with K° replaced by AKq — 
(JiÍÍj'q — Ï2K2Q) / (ji + Ï2) ■ Being the spectral contributions of the làser excess noise and phase diffusion concentrated 
in the low frequency region, it is worth noting that (see Eq. (§3-a)) 



AKq (0) = 1 - i-p^— ) — cos-0 



Eeff J 1-5* 



so that 



< (o) 
<4 (o) 



8x 
1-S 2 



1-S 2 



1 1 1 sin ^ • / 
cos ip cos ipp H sm %p p 

E eff 

cos ip sm yjp — — — cos tpj. 



E, 



eff 



(34) 



These expressions evidence the deleterious effects of crystal losses and detuning also in the difference spectra. While 
for a resonant device only the excess noise influences the spectrum, in a detuned device also the phase diffusion 
contributes to the noise in proximity of the origin. 

The plots of Fig. 8 show the difference spectrum Sd for an unbalanced resonant OPO (S — 0.05, -0 = 0) for 70 = 4 
(top) , 0.5 (bottom). The peak observed in the single beam spectrum for % = 0.5 (Figs. 6 and 7-bottom) is present 
also in the difference spectrum for the same damping ratio whereas for "f' = 4 (top) the difference spectrum keeps 
unchanged its sub-shot-noise character. 



IX. CONCLUSIONS 



The steady-state characteristics of the beams generated by an OPO depend on the decay constants of the cavity 
modes, the non linear crystal losses, the excitation parameter E, and the deviation from resonance (detuning). The 
last condition has been represented by an angle tp. While for resonant devices the beam amplitudes are proportional 
to \/ E — 1, in presence of detuning E has been replaced by an effective excitation parameter E . Similar ly, the phase 
difference between the driving field and the pump mode is a function of the detuning angle The effects of thcse 
changes propagate to the fluctuating parts of the twin beams amplitudes. 

Once linearized, the fluctuations of the twin beams are characterized by resonance freqüències which control the 

frequency responses K k \ 2 to the different noise sources (quantum noise of signal/idler/pump and phase-amplitude 
fluctuations of the pump làser). These functions have in general six poles fl r in the half plane Im(w) > 0. One, 
fio = 0, always lies at the origin, a second one, íii = iwx, is located on the positive imaginary axis, while the remaining 
four ones may be purely imaginary or complex; when complex (ÍÍ2 = — ^3, ^4 = — ^5) they are characterized by two 
damping constants W2,W4 and freqüències 1^2,^4- 

At resonance the K kl have only three poles, different from the those of K kl . Their locations in the complex plane 
depend on the excitation parameter E, the cavity damping coeficient 7 at the pump frequency, and the mismatch 
parameter S. For given 8 the couple of parameters (E, j' ) corresponding to purely imaginary resonances forms, on the 
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E—j' Q plane, a connected region delimited by the straight line E = 1 and a curve representing the maximum excitation 
£ maï for which the poles are all imaginary. For 7g comprised in a particular interval the poles are imaginary for 
1 < E < Ei and E 2 < E < EVna.y- Outside this region the two poles are complex conjugate and the spectrum exhibits 
a relaxation oscillation peak . 

The amplitude quadrature responses Kj k to the various noise terms are regular for lu — ► , contrariwise to the phase 

quadratures, KJ^ and K?y which present a pole. The difficultics of this pole have been bypassed by displacing it from 
the origin by a small quantity e, and letting e — > at the end of the calculations. So doing the complex fluctuation fj,—icf> 
of each mode emerges as proportional to the annihilation operator ò of a squeezed thermal oscillator (/i — i<j> = ha). 

By exploiting the Gaussian statistics of the fluctuations second and fourth order correlations have been expressed by 
suitable combinations of the frequency responses and their residues. So doing, it has been filled the gap between the 
evaluation of the variance and covariance of the fields carried out by Reid and Drummond p3| by means of classical 
stochastic equations, limited to resonant and symmetric OPO, and that based on quantum Langevin equations. The 
operator representation offers the opportunity of using the commutators of the field variables at different times for 
characterizing the respective correlations and gaining further insight in the detuning effects. In fact, the complex 



K°* M + IQ* (-w) 



2 



,0. 



2 



damping coeficients make asymmetric the frequency responses in prcscncc of detuning, 
thus implying that fj,j (t) and (f>j (i) at different times do not commute as does the instantaneous frequency with the 
phase at the same time (/ (f>j (0) <fij (0)^ ^= O^j . 

The detuning interplays with pump excess noise and phase diffusion by increasing the spectral broadening of the 
twin beams and enhancing the low frequency part of the spectra. These negative effects are mitigated by a well 
balanced configuration of the OPO. 

While triply resonant devices reduce the threshold levels by allowing the investigation of OPO quantum properties 
for large excitation parameters, their spectra differ notably from those typical of the adiabatic regime. The single 
beam spectral density starts, for lu = 0, from below the SNL, rapidly increases by reaching a peak in correspondence 
of the relaxation frequency, and then decays monotonically toward the SNL. The height of the peak increases notably 
as the pump decay constant decreases with respect to those of the signal/idler. These oscillations tend to obscure the 
cxpcctcd single bcam squcczing, by confming it to a low frequency region where the fluctuation spectrum is mostly 
conditioned by the làser excess noise and phase diffusion. 

X. APPENDIX A 

Ak (lo) and Bk (lu) depend on the functions Ak (lo) 



A ï A o ,__ A _, Aí , A* 



2 



Aj ~ Aj + A f A + l A*A* + 1' " U "" ü " r A a AÍ-l" r A 1 A*-l 

1 11 1 

B j — ~t — i - H 1 — i , Bo = + ! ï 

A,vAq + 1 A^At + l' A 2 A*-1 AiAl-1 

while the Ck (w) are linear combinations of the noise Fourier transforms Z' k (see Eqs. ([ïq)) 

- A^-Z> g? + A*2g A\Z' 2 + Z'Í Ajz[ + Z> 2 

C] Z Í + A.„A n + l + At. At , 1 ' 60 - Z ° A„AÍ_1 A~7t 



AJ,A5 + 1 A 2 Aí-l A^l 

For a detuned system (tp ^ 0) Ak differs from the adjoint A k , while Bk coincides with BÏ. 
Next, introducing the functions 

A k ± B k 



Fk± — 



AkAl-B* 

the frequency responses K k \ 2 defined by Eqs. @ 

are given by 



K°. = F* K°, = Fj - A °~ Fj - 

" 3 ~ ' " AjvAo + 1 

(| ___F}_A r + Fj _ „ F _+4_4 



K u — ■>— J K' 

K >°- A f A + l ' K °i- Aj-AÍ 

J J 
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and 



" ~ ■>+ ' »' " A,vA + 1 ' j0 ~ A,vA + 1 



Expressing Fk± and A& as functions of Cj yields 



K°- 

33 



K 



i + sa é* ^ + k^ü* + k°P& + + kT ] 



cos ip D' 

l-í^3 + ^(4)-2 + ^(3)- +ií 0( 2 ) 



n cos 2 ip D' 



D' 



and 



K 2 

33 



(l + Sj) *** ^ + gj°V + jffiV + gffV + KJ^Cj + 



cos ijj D' 
l-P^ + KJ^Cj* + K$ V + + 



JJ cos 2 V -D' 



Ki Q = %{E BÍf -l) (1 + 6. 



jO - 10 \^e}) 7 ^ ' 3 J 



ui 



with D' (w) defined in (|ï|) 



At resonance the -K"J 2 reduce to 



-,2 



-i l+fo-íj S-f B 1-í 



w + «70 [E - 2) 



^°o = -70(^-1) (l + íi) 



Cj — i2 ( 1 — 6, 



D 



+ 



and 



K?. =-i[l + S 



u 2 -i(l + %- ~8j) Q - %E (l - S s ) i 



^ = (l-<5 



D- Cj 



KJ Q = -% (E - 1) (1 + 5j) — 



In the adiabatic case (% ^> 1) 



A 30 



w - ££7 ( i - 
+ ^ . 

V ' Cj 2 -2iEJj-4(E-í)(l-S 2 



(E - 2) (l - ~S 2 ) - ^ 

V ' Cj 2 - 2íECj - 4 (E - 1) M - S 2 ) 



. . W — £2 ( 1 — S; 

(£7-1) (l + ^j- 



w 2 - 2iEuj - 4 (£7 - 1) 1 - ó 2 
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and 



K 2 - 

33 



33 



5, 



K 2 - 
-"70 ~ 



-H 1 

e(i-S 2 )^ 

\ ) LO — : 

-i(E-l)( 



lo — iE (í — Sj^J j 



lo — i2E LO 
1 

o — i2E üj 

1 



ü) - Í2E 



XI. APPENDIX B 



The coeficients of the characteristic polynomial are given by 



£)(5) 

£>(3) 
£)(2) 



= -i2(f + 2) 

= - (70 + 2) 2 - 4f - 45 2 tan 2 ^ - 4 (£ e// - 1) 7^ 
E eff % + 2 + 3 (E eff -!) + (! + E eff ) ~5 2 



*47Ó 
47Ó 



{(S e// - l) 2 + 3 (E eff - 1) + l}f + 2 (S e// - 1) - {(2 + f ) (£7 e// - 1) - y tan 2 ^} 



£)(i) = -i87^ e // (S e// - 1) (l - S 2 ) (37) 
with 70 and 5 respectively the normalized pump mode damping and mismatch factor 
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Figure captions 



1. Zeros fli, 1^2, of D' (u>) (see Eq. (l·q)), versus pump parameter E for ip = 0.3 and cavity parameters 7q 



2, S = 0.1. tui 



Im 



fil 



(dashed), cD 2 = Re 



and ui2 = Ini 



(dotted), and £4 = Re 



and 



u>4 = Im 



íï 4 



(continuous). Noticc that Re 



discussion illustrated in Fig. 2. 



Í2 2 ,4 



vanishes for E less a threshold value, in agreement with the 



2. Typical shape for 5 = 0.05 of the curve separating the points of the plane % — E, for which the three resonance 
freqüències of the amplitude quadratures /i are purely imaginary (below), from those with two complex conjugate 
roots (above). Relaxation oscillations occur only for points lying above the curve. 



3. Frequency Q2 = Re 



fi 5 



^Re Ò3 ^ (continuous) and damping coefficient 



W2 = Im 



2.3 



(dotted), of the 



relaxation oscillations at resonance (tp = 0), versus the pump parameter E (see Eq. fll9|)) for cavity decay rates 
% = 0.5, 2.5, 4.5,6.5 and 5 = 0. 



4. In-phase absolute squared transfer functions I (ü)) I , I K® 2 (lü) I ^ , and |ií° (ó))| z at resonance [tp = 0) vs 
normalizcd frequency for 7q = 1, § = 0.05, and pump parameters E — 2,3,4. 

5. Average phase commutator ([4>j (0) , <pj (r)]) versus r for % = 4, £ — 1.5, 4 and ip = 0.1. 

6. Single mode spectra for pump cavity damping 7q = 4 (top), 0.5 (bottom), 5 = 0.05, pump parameter E = 3, 
detuning phases tp = (sòlid), 0.25 (dashed), 0.5 (dot-dashed), and crystal losses X1/7 = 0.3, x n = xi/3. The 
spectra have been normalized to the shot-noise-level. 

7. Effects of the detuning on the transfer of pump excess noise to beam nuctuations. Plots parameters are as 
follow: 70 = 4 (top), 0.5 (bottom), S — 0.05, pump parameter E = 3, detuning phases ip = (sòlid), 0.25 
(dashed), 0.5 (dot-dashed), and crystal losses X1/7 = 0.3,>ío = ^i/3. 

8. Difference spectral density Sa for cavity damping 7Ó = 4 (top), 0.5 (bottom), 5 — 0.05, pump parameter E = 3, 
zero detunings (ip = 0), and crystal losses different for the pump and signal/idler ^1,2/7 = 0.3, xq = >fi/3. The 
spectra have been normalized to the shot-noise-level. 
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